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Abstract. We consider multipoint Pade approximation to Cauchy transforms of com- 
plex measures. We show that if the support of a measure is an analytic Jordan arc 
and if the measure itself is absolutely continuous with respect to the equilibrium 
distribution of that arc with Dini-smooth non-vanishing density, then the diagonal 
multipoint Pade approximants associated with appropriate interpolation schemes 
converge locally uniformly to the approximated Cauchy transform in the comple- 
ment of the arc. This asymptotic behavior of Pade approximants is deduced from 
the analysis of underlying non-Hermitian orthogonal polynomials, for which we use 
classical properties of Hankel and Toeplitz operators on smooth curves. A construc- 
tion of the appropriate interpolation schemes is explicit granted the parametrization 
of the arc. 
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1. Introduction 

Rational approximation to analytic functions of one complex variable is a most classical 
subject which has undergone many developments since C. Runge's proof that such an 
approximation is indeed possible, locally uniformly on the domain of holomorphy |47j . 
Let us quote for example the deep study of those open sets for which holomorphic func- 
tions that extend continuously up to the boundary can be approximated by rational 
functions on the closure |36} I60j , see comprehensive expositions and further references in 
the monograph |18j . In another connection, the achievable rate of convergence of rational 
approximants at regular points, when the degree goes large, also received considerable 
attention [HIl [HJ HHJ |32J US]- Meantime, rational approximation has become a tool in 
numerical analysis [3U1 [SS], as well as in the modelling and control of signals and 
systems, see for instance [2| [14 1 [8] [20 l [43] . 

From the constructive viewpoint , a great deal of interest has been directed to those 
rational functions of typer] (m, n) that interpolate a given function in m -I- n -I- 1 points. 



rational function is said to be of type (m, n) if it can be written as the ratio of a polynomial of degree 
at most m and a polynomial of degree at most n. 
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counting multiplicity. These are the so-called multipoint Pade approximants [6], which 
subsume the classical Pade approximants that interpolate the function in a single point 
with multiplicity m + n+1 [41j . Beyond the natural appeal of such an elementary linear 
procedure, reasons why Pade approximants have received continuing attention include 
their early number theoretic success when applied to certain entire functions |26ll50ll?T] . 
the remarkable behavior of diagonal (i.e. of (n,n)-type) Pade approximants to Markov 
functions [M! ^SS], the de Montessus de Ballore theorem and its generalizations on the 
convergence of (m, n)-type Pade approximants for fixed n |35] as well as the solution to 
the "disk problem" for diagonal sequences [22 that give rise to extrapolation schemes 
for analytic functions, the Nuttall-Pommerenke theorem on convergence in capacity of 
diagonal Pade approximants to functions having singular set of zero capacity |37[ I44j . 
the numerical use of Pade approximants for initial value problems and convergence ac- 
celeration |13l 1251 1271 112] , their connections to quantum mechanics and quantum field 
perturbation theory [71 (5^, and the research impetus generated by the so-called Pade con- 
jecture on the convergence of a subsequence of the diagonal sequence [S], in the largest 
disk of meromorphy, that was eventually settled in the negative |31j . The reader will 
find a detailed introduction to most of these topics, as well as further references, in the 
comprehensive monograph [6]. 

In the present paper, which deals with the convergence of multipoint Pade ap- 
proximants to functions defined as Cauchy integrals over a compact arc, the relevant 
approximants are the diagonal ones with one interpolation condition at infinity for they 
are those producing the (generically simple) zero assumed by the function at infinity. 
The first class of Cauchy integrals for which diagonal Pade approximants were proved 
convergent is the class of Markov functions, that is, Cauchy transforms of positive mea- 
sures on a real segment [5HI23]. In fact, under mild conditions, the diagonal multipoint 
Pade approximants to such functions at a conjugate-symmetric system of interpolation 
points converge strongly (i.e. locally uniformly in the domain of analyticity) , as the num- 
ber of interpolation conditions goes large. The working pin of this result is the intimate 
connection between Pade approximants and orthogonal polynomials: the denominator of 
the n-th multipoint Pade approximant is the n-th orthogonal polynomial with respect 
to the measure defining the Markov function, weighted by the inverse of the polynomial 
whose zeros are the interpolation points (this polynomial is identically 1 for classical Pade 
approximants). This is a key to sharp convergence rates, using the asymptotic theory of 
orthogonal polynomials with varying weights |23l I56j . 

This result has been generalized in several directions, in various attempts to develop 
a convergence theory of rational interpolants to more general Cauchy integrals. In particu- 
lar, the strong convergence of classical diagonal Pade approximants to Cauchy transforms 
of comp/ea;-?;a/Me(i functions on a segment was established in 10, 37, 39J, when the density 
of the function with respect to the arcsine distribution on the segment is smoothly invert- 
ible, and in |33j under the mere assumption that the function has continuous argument 
and vanishes almost nowhere. It is interesting that the first three references proceed via 
strong asymptotics for non-Hermitian orthogonal polynomials on a segment, whereas the 
last one relies on different, operator theoretic methods. Multipoint Pade approximants 
to Cauchy transforms of functions with non-vanishing analytic density with respect to 
the arcsine distribution on a real segment were in turn proved strongly convergent under 
quantitative assumptions on the near conjugate-symmetry of the interpolation scheme 
the proof rests on strong asymptotics for non-Hermitian orthogonal polynomials 
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with varying weight that dwell on the Riemann-Hilbert approach to classical orthogonal 
polynomials [15] as adapted to the segment in [28]. Let us mention that the occurrence 
of zeros in the density can harm the convergence, because some poles of the approximant 
may no longer cluster to the segment in this case, but if the density has an argument of 
bounded variation the number of these "spurious" poles remains bounded and one can 
still get convergence in capacity even for substantial zeroing of the density ^ . 

The case of Cauchy integrals over more general arcs than a segment turns out 
to be harder. In the series of pathbreaking papers [SU [531 IMl classical diagonal 
Fade approximants to functions with branchpoints were shown to converge in capacity 
on the complement of the system of arcs of minimal logarithmic capacity outside of 
which the function is analytic and single- valued. This extremal system of arcs, often 
called nowadays a symmetric contour or an S-contour, is characterized by a symmetry 
property of the (two-sided) normal derivatives of its equilibrium potential, and the above- 
mentioned convergence ultimately depends on a deep potential-theoretic analysis of the 
zeros of non-Hermitian orthogonal polynomials over S'-contours. Shortly after, the result 
was extended to multipoint Fade approximants for Cauchy integrals of continuous quasi- 
everywhere non-vanishing functions over S'-contours minimizing some weighted capacity, 
provided that the interpolation points asymptotically distribute like a measure whose 
potential is the logarithm of the weight [24. With these works it became transparent 
that the appropriate Cauchy integrals for Fade approximation must be taken over S'- 
contours, and that the interpolation points should distribute according to the weight that 
defines the symmetry property. Subsequently, strong convergence of diagonal multipoint 
Fade approximants was obtained in this setting for analytic S'-contours and non- vanishing 
analytic densities with respect to either the equilibrium distribution of the contour or the 
arclenth, under certain quantitative assumptions on the convergence of the interpolation 
points 3 . Earlier results in a specific case can be found in |57j . 

Surprisingly perhaps, the natural inverse problem whether given a system of arcs, 
say S, there exists a logarithmic potential whose exponential defines a suitable weight 
making S into an S-contour, and whether the interpolation points can be chosen ac- 
cordingly to produce strong convergence of multipoint Fade approximants to Cauchy 
integrals on S, was apparently not considered. The goal of the present paper is to fill this 
gap when S is smooth. More precisely, we prove that a smooth arc is an S'-contour for a 
weight whose logarithm is the potential of a positive measure supported disjointly from 
the arc if, and only if the arc is analytic (c/. Theorem [T|. When this is the case, we further 
show there exists a scheme of interpolation points producing strong convergence of the 
sub-diagonal multipoint Fade approximants to Cauchy integrals on the arc, provided the 
density of the integral with respect to some positive power of the equilibrium measure of 
the arc is Dini-smooth and non- vanishing (c/. Theorems |2] and |4]) . This result seems first 
to ascertain convergence of an interpolation scheme for Cauchy integrals of fairly gen- 
eral densities (note analyticity is not required) over a rather large class of arcs (namely 
analytic ones). Finally, we show that the non- vanishing requirement on the density can 
be relaxed a little (c/. Theorems |3] and |4]) , thereby providing an initial example, over 
an arbitrary analytic arc, of a family of functions whose zeroing at isolated places does 
not destroy the strong convergence of multipoint Fade approximants to their Cauchy 
integral. 

Our proofs differ somewhat from those commonly encountered in the field. In par- 
ticular, we translate the symmetry property, which is of geometric character, into an 
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analytic one, namely the existence of a sequence of pseudo-rational functions tending to 
zero off the arc and remaining bounded on the latter. Choosing the interpolation points 
to be the zeros of these pseudo-rational functions, we transform the non-Hermitian or- 
thogonality equation satisfied by the denominator of the multipoint Pade approximant 
into an integral equation involving Hankel and Toeplitz operator on the arc S. Then, the 
bounded behavior of our pseudo-rational functions teams up with classical compactness 
properties of Hankel operators to produce strong asymptotics for the non-Hermitian or- 
thogonal polynomials under study. This in turn provides us with locally uniform error 
rates for the approximant on the complement of S. From the technical point of view the 
paper can be seen as deriving strong asymptotics for non-Hermitian orthogonal poly- 
nomials with complex densities with respect to the equilibrium distribution of the arc 
under minimum smoothness assumptions as compared to those currently used in the 
Riemann-Hilbert methods |35l E] . 

The paper is organized as follows. In the next section we formulate our main results 
(Theorems [l]-[4]). Section |3] contains the necessary background material for the proofs 
adduced in Section|4] The last section contains an example of a class of contours for which 
the behavior of orthogonal polynomials is "highly similar" to the case when F = [—1,1]. 



2. Main Results 

Let be a rectifiable Jordan arc with endpoints ±1 oriented from —1 to 1. Set 

w{z) = w(F; z) := \/ z^ — 1, 'w{z)/z — > 1 as z oo, (2-1) 

which is a holomorphic function outside of F. Then w has continuous traces (boundary 
values) from each side of F , denoted by w+ and w~' (e.g. w'^ is the trace of w taken 
from the left side of -F as -F is traversed in the positive direction). In this paper we 
consider polynomials g„, deg(q„) < n, satisfying non-Hermitian orthogonality relations 
with varying weights of the form 

t^qn{t)Wn{t)^^^Q, j-0,...,n-l, (2.2) 

together with their associated functions of the second kind, i.e. 

i?„(.):=l / ^J^^J^ zeD:^C\F, (2.3) 
m Jp t — z w+{tj 

where C is the extended complex plane and {w„} is a sequence of complex- valued func- 
tions on F. Our main goal is to show that the polynomials g„ possess so-called strong 
(Szego-type) asymptotics under the assumption that Wn — hn/vn, where is a com- 
pact family of sufficiently smooth functions on F and t;„ are polynomials of respective 
degrees at most 2n, subject to certain restrictions. The desired result can be expressed 
as 

(2/^)"(g„S^J = l-f o(l), ^■.= z + w, (2.4) 
locally uniformly in D, where Sw„ is the Szego function of Wn (see Section . As a 



consequence of (2.4 1, we establish the uniform convergence of multipoint diagonal Pade 
approximants to Cauchy transforms of complex- valued measures of the form hdt/w^, 
where /i is a sufficiently smooth function on F and the Pade approximants interpolate 
the corresponding Cauchy transform at the zeros of the polynomials f„. 
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2.1. Symmetry Property 

Let F be as above and D be its complement in C. In what follows, we always assume 
that the endpoints of F are ±1. Define 

if{z) := z + w{z), z e D, (2.5) 



where w is given by (2.11. Then is a non- vanishing univalent holomorphic function in 



D except for a simple pole at infinity. It can be easily checked that 

— —w^ and therefore (p'^'p^ = 1 on _F. (2-6) 

Let e E D. Define 

r{e; z) := ^^^"^ , |e| < oo, and r(cx); z) z <E D. 

Clearly, r(e; •) is a holomorphic function in D with a simple zero at e and non-vanishing 
otherwise. Moreover, it follows from (2.6) that unrestricted boundary values r^(e; •) exist 



everywhere on F and satisfy 

r+(e;i)r-(e;i) = 1, t e F. (2.7) 

Let now S' :— {En}neN be a triangular scheme of points in D, i.e. each En consists 
of 2n not necessarily distinct nor finite points contained in D. We define the support of 
S as supp((?) := C\n&i^k>nEk C C. Clearly, the support of any weak* limit point of 
the counting measures of points in En is a subset of supp(fS'). Hereafter, the counting 
measure of a finite set is a probability measure that has equal mass at each point counting 
multiplicities and the weak* topology is understood with respect to the duality between 
complex measures and continuous functions with compact support in C. To each set En 
we associate a function r„ by the rule 



rn{z) Jl r{e-z), z € D. (2.8) 

Then r„ is a holomorphic function in D with 2n zeros there and whose boundary values 
on F satisfy r+r~ = 1. Among all possible schemes S', we are interested only in those 
that induce the following symmetry property on F. 

Definition (Symmetry w.r.t. S"). Let F be a rectifiable Jordan arc and S be a triangular 
scheme of points in D. We say that F is symmetric with respect to S' if the associated to 
S functions rn satisfy |r^| = 0(1) uniformly on F and Tn = o(l) locally uniformly in D. 

Note that the boundedness of |r^| from above implies their the boundedness away 



from zero by (2.7) 



As we shall see, this definition is, in fact, closely related to the classical definition 
of the symmetry of an arc in a field (more generally, of a system of arcs) . The latter is 
based on a number of potential-theoretic notions all of which can be found in |49) . 



Definition (Symmetry in a Field). Let F be a rectifiable Jordan arc and let q be a har- 
monic function in some neighborhood of F. It is said that F is symmetric in the field q 
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if the following partial derivatives exis^ and are equal: 

d{U^ + q) d{U^ + q) 

Qn+ = g„- supp(A), (2.9) 

where A is the weighted equilibrium distribution in the field q, is the logarithmic 
potential of X, and ri^ are the one-sided unit normals to F. 

The symmetry property turned out to be vital in the investigation of the rates of 
best approximation of functions with branch points by rational interpolants (multipoint 
Pade approximants) |53l 155] . Given a continuum E with connected complement and a 
function / holomorphic in a neighborhood of E, analytically continuable except over a 
compact set Eg, cap(Fo) = 0, where cap(-) is the logarithmic capacity, then there exists 
an extremal (maximal) domain D such that the condenser (E,E), E := dD 3 Eq, has 
minimal condenser capacity among all domains in which / has single- valued holomorphic 
continuation. Moreover, if the continuation of / over C\Eq is not single-valued then the 
boundary E of the extremal domain is a union of a system of open analytic arcs and 
a set of capacity zero, and it can be characterized as the unique contour containing Eq 



and satisfying (2.9 1 with q = — JJ"^, where i^e — vp is the equilibrium distribution for 
the minimal energy problem for signed measures of the form ge — cffi where ap and 
ap range over all probability Borel measures on E and E, respectively. The best rate 
(in the n-th root sense) of approximation to f on E is achieved by rational interpolants 
corresponding to a triangular scheme whose points are asymptotically distributed as ve- 
The counting measures of the poles of such interpolants weakly converge to vp and the 
interpolants themselves converge to / in capacity in D. 

Dwelling on the work of H. Stahl discussed above, A. A. Gonchar and E. A. 
Rakhmanov |24j extended the definition of a symmetric contour to general harmonic 
fields in order to give a complete formal proof of the "Magnus conjecture" on the rate of 
approximation of the exponential function on the positive semi-axis[^ In [24^ it is assumed 
that E and q are such that supp(A) is a tame set, i.e. the intersection of a sufficiently small 
neighborhood of quasi every point (q.e. means up to a set of zero capacity) of supp(A) 
with the support itself is an analytic arc. It is easy to check that in this case partial 



derivatives in (2.9 1 do exist. It was also observed that in the setting of rational interpola- 
tion to functions as above, one can take interpolation points asymptotically distributed 
like any Borel measure v on E. This will define a unique contour E D Eq that satisfies 



(2.9) with q — —If^ and the interpolants will converge to the approximated function in 
capacity in C \ while the counting measures of their poles will weakly converge to i/, 
where p is the balayage of v onto E. 

Since our main interest lies with rational interpolation, we shall be concerned only 



with arcs satisfying (2.9 1, where q = —V^ for some Borel measure v with compact 
support. In this case much milder assumptions on E are sufficient in order to have well- 
defined partial derivatives (see the statement of the next theorem) . Necessarily, for such 
fields, the symmetric arcs turn out to be analytic as apparent from the following theorem. 



■^Normal derivatives are understood in the strong sense, namely, if the tangent to F exists at t and are 
the unit normals from each side of F then the limits of the Euclidean scalar products (V((7'*' +q){y),rV^) 
exist as y approaches t along , respectively. 

^This conjecture was formerly known as "1/9 conjecture", which was shown to be false by H.-U. Opitz 
and K. Scherer in |40| . The correct constant was identified by A. P. Magnus in |32l . however, his proof 
was not entirely formal. 
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Theorem 1. Let F be a rectifiable Jordan arc such that for x — ±1 and all t ^ F 
sufficiently close to x it holds that \Ft^x\ < const. |a; — t\^ , (3 > 1/2, where \Ft,x\ is the 
arclenth of the subsarc of F joining t and x and const, is an absolute constant. Then the 
following are equivalent: 

(a) there exists a triangular scheme of points S , supp((o') C F) , such that F is symmetric 
with respect to S" ; 

(b) there exists a positive compactly supported Borel measure v , supp(i^) C F , such that 
F is symmetric in the field ; 

(c) F is an analytic Jordan arc, i.e. there exists a univalent function p holomorphic in 
some neighborhood of [—1, 1] such that = 1, 1]); 

The above theorem covers only the case where supp((f ) is disjoint from F. The 
authors do not know whether non-analytic arcs can be symmetric with respect to a 
triangular scheme when the support of the latter does intersect the arc. 

We point out that the proof of Theorem [T] is constructive. In other words, for a 
given analytic arc F, a suitable scheme ca can (in a non-unique manner) be explicitly 
written in terms of the function p that analytically parametrizes F. Each such scheme 
(o gives rise to a suitable measure v simply by taking the weak* limit of the counting 
measures of points of the sets in S". 

2.2. Strong Asymptotics for Non- Vanishing Densities 

Let K he a, compact set in C and denote by C{K) the space of continuous functions 
on K endowed with the usual supremum norm || • \\k. For h G C{K), set uj^ to be the 
modulus of continuity of h, i.e. 

LOhir) -.^ max \h{ti) - h{t2)\, t e [0,dia,ia{K)], 

\tl-t2\<T 

where diam(i^) := maxf^ \ti — t2\. It is said that h is Dini- continuous on K if 

f ^hir) 

/ dr < oo. 

J[0,diam(if)] 

We denote by DC{K) C C{K) the set of Dini-continuous functions on K and by DCuj{K) 
the subset of DC{K) such that w?i < w for every h € DC{K), where oj is the modulus of 
continuity of some function in DC{K). In what follows, we shall employ the symbol "*" 
to indicate the nowhere vanishing subset of a functional set (for instance, C*{K) stands 
for the non- vanishing continuous functions on K). 

Let h G DC* {F) and log h be an arbitrary but fixed continuous branch of the 



logarithm of h. Then it is easily verified (see Section 3.3 1 that the geometric mean of h, 
i.e. 

G.:=exp{/^logMOj^}, (2.10) 

is independent of the actual choice of the branch of the logarithm and is non-zero. More- 
over, the Szego function of h, i.e. 

^ / N ( wiz) f logh(t) idt 1 f ^ , . idt '] ^ ^ 

S^{z) :^ cxp ^ / ^ \' - - / logfeft^— , zeD, 2.11 

1^ 2 Jp z-t TTW+{t) 2 Jp nw+{tj ) 

is also independent of the choice of the branch (as long as the same branch is used in 
both integrals). In fact, it is the unique non- vanishing holomorphic function in D that 
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has continuous unrestricted boundary values on F from each side and satisfies 

h = GhS^S^ on F and Shioo) = 1. (2.12) 

To state our next theorem we need one more notion. Let X be a Banach space 
and y be a subset of X. We say that a family {h„} C K is relatively compact in Y if 
any sequence from this family contains a norm-convergent subsequence and all the limit 
points of {hn} belong to Y. 

With the use of the previous notation, we formulate our first result on strong asymp- 
totics of non-Hermitian orthogonal polynomials. 

Theorem 2. Let F be an analytic Jordan arc connecting ±1 that is symmetric with respect 
to a triangular scheme of points S' = {En}n€n- Further, let {qn}n£N be a sequence of 
polynomials satisfying 

t^qnit)wn{t) = 0, j = 0,...,n-l, 

with Wn = hn/vn, whcrc {hn\ is a relatively compact family in DC*{F) and u„ are 
m.onic polynomials with zeros at the finite points of En . Then, for all n large enough, 
polynomials qn have exact degree n and therefore can be normalized to be monic. Under 
such a normalization, we have that 

Rnw Z [lto(l)hS„ locally umformly zn D, (2.13) 



where Sn '■— i'^/v)"'Sw„, 7n 2 Gw^, and _R„ was defined in (2.3 I 



Qn = (l + d„)/5+ + (l + d+)/5„ , . 



\dt\ —tQasn 



where d^ € C{F) and satisf'^ 



F 



for any p G [1, oo). Furthermore, the following limit takes place 

^M^Wrfi A (2.15) 

-fnW+{t) W+{t)' ^ ' 

where dt is the differential along F and "-^ " stands for the weak* convergence of mea- 
sures. 

The theorem is stated for analytic arcs only because in the view of Theorem [T] 
this will be the case if the scheme S' does not touch F, which is a standard setting in 
interpolation. In fact, a careful examination of the proof of Theorem |2] shows that F 
could simply be a Carleson (Ahlfors-regular) curve. However, the authors do not have 
any convincing evidence suggesting that such arcs could be symmetric with respect to 
some triangular scheme unless they are analytic. 



*From classical estimates on singular integrals 1171 the functions (ij„ui)^ are, in fact, continuous func- 
tions on F; this also can be seen, for instance, from equation ( |3.5[ l below. However, | |2.14[ l does not 
contain explicit information on the pointwise behavior of (Rnw)^ as the smallness of dn is claimed in 
LP norm only and pays no particular attention to the values of d„ on a set of zero linear measure of F. 
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2.3. Strong Asymptotics for Densities with Some Vanishing 

In this section we consider an extension of Theorem [2] to densities that may vanish. As 
mentioned in the introduction, such an extention is important since it provides the first 
example of strong asymptotics with vanishing density. The method used in the proof is 
essentially the same and is presented separately solely for the clarity of the exposition. 

The purpose of the following construction is to introduce a factor having zeros in 
the density of the measure and whose vanishing is sufficiently weak so that the proof of 
Theorem |2] still goes through with minor modifications. 

Denote by Fed the closed subarc of F with endpoints c and d. As F is compact 
there exists y e (— oo,— 1] such that F n (— oo,y) = 0. Chose Faux to be any smooth 
Jordan arc, Faux HF = {— 1}, that links y and —1. Let x E F, define 6^ to be the branch 
of arg(- — a;) that is continuous in C\ {(— oo, y] Ui^aux UF_i_2;} and tends to zero when the 
variable approaches infinity along the positive real axis. For definiteness, we set 9x{x) to 
be the limit of 9^ along F^^i- Then 6^ has continuous traces 0^ on F-i^^ oriented from 
— 1 to X and 0^{t) = 9~{t) + 2tt, t E F^i ^ \ {x}- Clearly, the functions 6^ and O^^p^ ^ do 
not depend on the choice of Faux ■ 

For any a E (0, 1/2] and a; G F set 

;.(a,:.;t):=|2(t-x)p4 '"Pg*"?£)/' l"" s x (2-16) 

^ ' ' ^ I ^ ^1 \ exp{2ia(6'+(<) -tt)}, tEF_i,^\{x}. ^ ' 

Then h{l,x;t) — 4{t — x)'^, t E F, the argument of h{a,x; ■) is continuous on F, and 
n{a,x;t) = |2(t-x)|2" when F= [-1,1]. 

Let hhe a, function on F for which Sh is well-defined. Set / . It is shown 

that functions are continuous on F and satisfy c^(±l) = 1 whenever 



3.3 



in Section 

h E DC*{F). Moreover, it is also shown that the functions c^j-q, ^..-| are continuous on 



F \ {x} and have moduli t E F \ {x}, where (p was defined in (2.51. They 

"h{a,x;-)i 



possess the one-sided limits along F at x, denoted by cJ, . Jx*), and 



4ia..x,)i^'') - 4(a,.,)(^') = 2sin(a7r) \ip^{x)\'" , (2.17) 

where we make the convention that 

c+ ,(l+):=c^, ,(1~) and c~, ,(-1^) c+, ^(-l+). 

Now, let Fq C F be a finite set of distinct points. We associate to each a; G Fq some 
ax E (0, 1/2) and define 

h{t) = h{Fo;t) -.^ Yl h{a,,x;t) Y[ h{aj2,x;t), (2.18) 

xGFo\{±l} a:eFon{±l} 

t E F. Then the following theorem takes place. 

Theorem 3. Assume that 

• F is an analytic Jordan arc connecting ±1 that is symmetric with respect to a 
triangular scheme of points S = {F„}„gN and let w„ be the monic polynomial with 
zeros at the finite points of En; 

• the functions r„, associated to S via ( |2.8| , are such that {|r„ o iy9~^|}„gN is a rela- 
tively compact family in DC(T), where T is the boundary ofip{D); 

• {^n}neN is a relatively compact family in DC^{F); 
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• Fq C F is a finite set of distinct points, Ux G (0,1/2), x G Fq, and h, given by 
(2.18 1, is such that 

limsup|(cR/i„r„)±(a:+) - {cnh^r,,)^{x')\ < 2\\Q\\-\ (2.19) 

n — >oo 

X G Fq, where 1 < ||Q|| is the norm of the outer Cauchy projection operator on 
L'^{T) defined m (|312|. 



If {qn}neN sequence of polynomials satisfying (2.2 1 with = hhn/vn, then the 
polynomials have exact degree n for all n large enough and therefore can he normalized 
to he monic. Under such a normalization (2.131 and (2.141 hold with dn satisfying 

\dnit)\' 



lim sup 



\dt\ < oo. 



Let us make several remarks. As in the case of Theorem |2] the requirement of 
analyticity of F can be weakened. It is enough to assume F to be piecewise Dini-smooth 
without cusps. This condition is sufficient for Section |3.6| to remain vahd, which is the 
main additional ingredient in the proof of Theorem [3] as compared to the proof of Theorem 

12 

As clear from the definition of Q, the lower bound in the inequality ||Q|| > 1 is 
achieved when F = [—1, 1]. 

To obtain an explicit upper bound for the numbers ax, x G Fq, for which the 
theorem holds, rewrite ( |2.19 1 in view of (2.171 as 



sin(aa;7r) k5'''(x) 



limsup|(cft„r„)±(a::)| < ||Q|1 



Thus, since 2ay <1, y € Fq, (2.191 is satisfied for example when 

sin(a:,7r) < ||Q||-i min{|(^-(x)|^, |(/3+(x)|^}liminf |(c,,„r„)±(a;)| 

n — >oo 

where L is the number of elements in Fq. In particular, if = [— 1, 1], the functions /i„ 
are positive, and the sets E„ are conjugate-symmetric, then \ip^ \ = |cj | = |r^| = 1 and 
(2.191 simply reduces to the initial condition 2ax < I, x E Fq. 

2.4. Multipoint Pade Approximation 

Let /i be a complex Borel measure with compact support. We define the Cauchy transform 
of /X as 

d^i{t) 



z-t 



z G C \ supp(/i). 



(2.20) 



Clearly, is a holomorphic function in C \ supp(/i) that vanishes at infinity. 

Classically, diagonal (multipoint) Pade approximants to are rational functions 
of type (n, n) that interpolate at a prescribed system of 2n+ 1 points. However, when 
the approximated function is of the from (2.201, it is customary to place at least one 
interpolation point at infinity. More precisely, let S = {£"„} be a triangular scheme of 
points in C \ supp(/z) and let w„ be the monic polynomial with zeros at the finite points 

of £;„. 



Definition (Multipoint Pade Approximant). Given of type (2.201 and a triangular 
scheme S , the n-th diagonal Pade approximant to associated with S is the unique 
rational function n„ = Pn/q-n satisfying: 
• degp„ < n, degg„ < n, and (j„ ^ 0; 
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• {(ln{z) f f_i{z) - Pn{z)) /v„{z) IS analytic inC\ supp(^); 

• iQn{z)fi_^{z) - Pn{z)) /Vn{z) = O (l/z"+l) flS Z -> OO. 

A multipoint Pade approxiniant always exists since the conditions for p„ and qn 
amount to solving a system of 2n + 1 homogeneous linear equations with 2n + 2 unknown 
coefficients, no solution of which can be such that = (we may thus assume that q„ is 
monic); note that the required interpolation at infinity is entailed by the last condition 
and therefore n„ is, in fact, of type (n — l,n). 

The following theorem is an easy consequence of Theorems [l] and [2] 

Theorem 4. Let F be an analytic Jordan arc connecting ±1. Then there always exist 
triangular schemes such that F is symmetric with respect to them. Let S be any such 



scheme and be given by (2.201 with 



idt 

dMi)=A(i)— TTTT, ii&DC^F), supp(/i)=F. (2.21) 
Trw^[t) 

Then the sequence of diagonal Pade approximants to associated with S , {!!„}, is such 
that 

(ffj — n„)i(; = [2Gf^ + o(l)]S'?r„ locally uniformly in D. (2.22) 

The proof of Theorem [4] combined with Theorem [3] yields the following. 

Corollary 5. Let F an analytic arc connecting ±1 that is symmetric with respect to a 
triangular scheme S'. Further, let h € DC*{F), Fq <Z F be a finite set of distinct points, 



and h{F(); ■) be given by (2.181 for some choice of the parameters < 1/2, x ^ Fq. Lf 



functions r„, associated to S , are such that {|r„ o ip ^\} is relatively compact in DC{T), 



r dif{D), and (2.191 is satisfied with hn = h and h — h(Fo; ■), then the conclusion of 



Theorem^ remains valid for f^ with (i = hh. 

As a final remark of this section, let us mention that it is possible to consider in 
Theorem |4] the Radon-Nikodym derivative of /i with respect to the equilibrium distribu- 
tion on F. The later will be different from /i by a factor g{ip)ip, where <? is a non- vanishing 



function analytic across T. However, as apparent from (2.221, the adopted decomposition 
of /i is more convenient. 



3. Auxiliary Material 

In this section we provide the necessary material to proceed with the proofs of the main 



results. Section 3.1 describes properties of F, the boundary of (p{D). In Sections 3.2 and 



3.3 we constructively define Szego functions, show that formula ( |2.12 1 indeed takes place. 



and compute several examples that we use throughout the proofs. The main source for 
these sections is the excellent monograph [T7]. Sections 3.4 and 3.5 are concerned with 
the so-called Smirnov classes on Jordan curves and their content can be found in much 
greater detail in [11]. Finally, Section 3.6 carries onto analytic Jordan domains certain 
properties of inner-outer factorization of analytic functions on the unit disk. 
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3.1. Joukovski Transformation 

Let F be a rectifiablc Jordan arc connecting ±1, oriented from —1 towards 1, and let 
tf be defined by (2.51. Then has continuous injective extensions tp^ onto F that are 
conformal on \ {±1}. Denote 

r:=r+ur-, t^:=lp^{f). 



By what precedes, \ {±1} are open Jordan arcs with endpoints ±1 and by (2.6) 

r- = (r+)-i :={z: i/zer+}. (3.1) 

Observe that is inverse to J{z) :— (z + l/z)/2, the usual Joukovski transformation, in 
D = C \ F. Since J maps z and 1/z into the same point and (p^{z) = ±1 if and only if 
z — ±1, r+ n = { — l, 1} and therefore F is a Jordan curve. 

Assume, in addition, that for x = ±1 and all t S F sufficiently close to x it holds 
that \Ft^j;\ < const. |a; — t\^, (3 > 1/2, where |-Ft is the arclenth of the subsarc of F 
joining t and x and const, is some absolute constant. We shall show that this implies 
rectifiability of F. As (p^ are conformal on F \ {±1}, it is enough to prove that F has 
finite length around ±1. By (3.1 ), it is, in fact, sufficient to consider only F"*". Let s be a 
parametrization of F around 1 such that s(0) = 1 and for all 6 E [0,do], s{So) ^ —1, it 
holds that 

m := = / \s'{y)\dy < const.|l - s{5)f. 

As if'^ o s is a parametrization of F+ around 1, we get that 

\^+iF,iS),i)\ = I \{^^{s{y)))'\dy^ I "P"}^ \s'{y)\dy 
J[o,5] J[a,5] w+{s{y)) 

< const. / ' ^^'^ dy 

since (if^)' = ip^ /w^ , is bounded above on F, and |1 + s{y)\ is bounded away from 
zero on [0,(5o]. As f{y) = |s'(y)| and |1 — s(y)| > const.£^^'^{y) on [0, (5o], we deduce that 

< const. ^iT^'^y = consU'''/^^{S) < oo 

since 1 — 1/2/3 > 0. Clearly, analogous bound holds around —1. Hence, F+ and therefore 
F are rectifiable. 

Now, we show that F is an analytic Jordan curve whenever F is an analytic Jor- 
dan arc. In other words, we show that in this case there exists a holomorphic univalent 
function in some neighborhood of the unit circle that maps T onto F. Let U he a, neigh- 
borhood of T such that J{U) lies in the domain of p, where p is a holomorphic univalent 
parametrization oi F, F = 1, 1]). Define 



^z) 



i/p{p{j{z))), zeu+:=un: 

^(p(J(z))), zEU-:^U\l 
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Then $ is a sectionally holomorphic function on U \ T. Denote by the traces of $ 
from C/± on T. Then for r e T, we have 



l/^+(p(J(r))), Im(T)>0 
1/^-(p(J(t))), Im(T)<0 

^-(p(J(t))), Im(T)>0 
ip+{p{J{T))), Ini(r) <0 ~ ^ >■ 



Thus, $ is a holomorphic injective function on U that maps T into T and therefore T is 
an analytic Jordan curve. 

3.2. The Cauchy Integral Operator 

Let r be an analytic Jordan curve and denote by and D~ the bounded and unbounded 
components of the complement of F, respectively. Further, let (j) be an integrablc function 
on F. The Cauchy integral operator on F is defined as 

C0(z):=-^ / ^dr, z^T. (3.2) 

Clearly, Ccf) is a sectionally holomorphic function on C \ F and 

t e F, 



' \ A4t,z), z& D- 



where 

27rj Jp T — z 

It is an easy modification of [17, Sec. 4.1] to see that A^{t,z) is a continuous function 
of z for each t E T when (f> E DC{T). Therefore, in this case, C(j)^o± have continuous 
unrestricted boundary values (C0)* such that 

(C0)+ - (C0)- = (3.3) 



Moreover, C4> is the unique sectionally holomorphic function in C \ F that satisfies (3.3 1 
and vanishes at oo. 

Let now {(j)n} C DC^{T) be a uniformly convergent sequence. Denote by (p the limit 
function of and observe that ^ necessarily belongs to DC^{T). We want to show 

that {C(f>n)^ converge uniformly on T to {C<f>)^ , respectively. Since C is additive, we may 
suppose without loss of generality that (j) = 0. Then 

\<Pn{r)-Mt)\ \dT\ 
\T-t\ 27r 



|(C^„)-(i)| = \A^At,t)\< 



< h I "^ds^kj f + f l^d. 



[0,so] ^ V"'!"^*"] "'[«":«o], 

[0,s„ 



< ki I ^^^ds + 2fci(so - Sn)- 

"'fO,s„,l S 



where sq diam(F) and fci is a positive constant that depends only on F. Now, by 
choosing s„ in such a manner that s„ and |j0n||r/5n both converge to zero as n — s- cx), we 
obtain that (C0„)~ converges to zero uniformly on F. The convergence of (C(/>„)^ then 
follows from (3.3 1. 



14 



L. Baratchart and M. Yattselev 



Finally, suppose that T is such that r and 1 /t belong to T simultaneously. Assume 
in addition that G DC{T) is such that 0(r) = 0(l/r). Then 

C?!)(l/z) =C0(O) -C(/)(z), z^r, 

(C(^)± e C(r), and by ^\ 

(C0) + (r) + (C(/.) + (l/T) = 0(T)+C0(O), rer. (3.4) 

Suppose further that ±1 G F and denote F = «/(F). Then F is an analytic Jordan arc with 
endpoints ±1. Let ip G DC{F) and (/) = V ° >/• Clearly, G i:)C(F) and 0(t) = ?!)(1/t), 
T G F. Set 

As 

dJ(T) _ r -dr/r, T G F" 



w+(J(r)) \ dr/T, TGF+ ' 

where F* have the same meaning as in Section 3.1 and since w{(^) = (1 — z^)/(2z), 
C = J{z), z G , we have for t — J{f) that 



where we used the symmetry of 0, i.e. (/)(r) = (/)(1/t). Then we get from (3.5 1 and ( |3.4[ | 
that 

G C(F), i?+(±l) = i?-(±l), and R+ + R-^i; on F. (3.6) 
3.3. Szego Functions 

Let F be an analytic oriented Jordan arc connecting —1 and 1 and h G DC*{F). Since 
|Log(l + z)\ < 2\z\, \z\ < 1/2, where Log is the principal branch of the logarithm, we 
have that 

|log/j(<i) -log/i(t2)| = \Log{h{h)/h{t2))\ < 25^^^il^l^M, (3.7) 

mmp h 

whenever 2ujh{\ti — t2\) < mmp h. Thus, log/i G DC{F) for any continuous branch of 
the logarithm of h. Let now F be the preimage of F under the Joukovski transformation, 
J, and F* ~ ip^{F). Further, let (f) := log ho J. Then, changing variables as in (3.5), we 
have that 

C. ^ cxp {^logMO^} - cxp{X ^|.} ^ »p(C0(O)} ^ 

and does not depend on the branch of the logarithm we chose. Indeed, different continuous 
branches of the logarithm of h give rise to the same (j) up to an additive integer multiple 
of 27ri, and since C((/)+ c)(0) = C0(O) + c for any constant c, the claim follows. Further, 
(3.5) yields that 

5„(C)=exp{C(/.(2)-C</)(0)}, C^J{z)eD, zeD+. 

Hence, 5"^ is a well-defined non- vanishing holomorphic function in D that has value 1 at 
infinity and satisfies 

St{t)S^{t)^exp{<p{T)-Cm}^h{t)/Gf,, t^J{r)eF, (3.8) 
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by (|3^ and (|3^. Moreover, (|3^ tells us that 

S^eCiF) and S+ (±1) ^ (±1). (3.9) 

Now, if /ii,/i2 G DC*{F) then S'^i/ia = Sh^S^^ since there always exists fc e Z such that 

log(/iift,2) = log hi + log /i2 + 27rifc on 

for any continuous branches of the logarithms of hi and /i2. We shall refer to this property 
as to the multiplicativity property of Szego functions. It remains only to say that can 
be characterized as the unique holomorphic non-vanishing function in D that assumes 



the value 1 at infinity, has bounded boundary values on F, and satisfies (3.8 1. Indeed, 



this follows from the uniqueness of C(j> as the sectionally holomorphic function satisfying 



(3.3 1, vanishing at infinity, and having bounded boundary values. 

Employing the characterization described above, we can provide an explicit expres- 
sion for the Szego function of a polynomial. Namely, let En, be a set of 2n not necessarily 
distinct nor finite points in £), Vn be a polynomial with zeros at the finite points of £"„, 



and r„ be the function associated to En via (2.8 1. Then we get from (2.6) and (2.7 1 that 



^l-7^r^- (3-10) 



1 

Proceeding now with the investigation of the properties of Szego functions, we show 
that {S^ } form relatively compact families in C*{F) whenever is a relatively com- 
pact family in DC*{F). Without loss of generality we may assume that {/i„} converges 
uniformly on F to some non-vanishing function h. Fix an arbitrary branch of the loga- 
rithm of h. Clearly, there exists a choice of the branch of the logarithm of /i„ for each n 
such that the sequence {log/i„} converges uniformly to log/i. Moreover, it follows from 



(3.7 1 that {log/i„} C DCi^{F). Here, we slightly abuse the notation by using the same 
symbol uj despite the fact that the modulus of continuity that majorizes wiog/i^ is dif- 
ferent from the one that majorizes lou^ . As only the existence of such a majorant is 
significant to our considerations, this should cause no ambiguity. Thus, C DC^iT), 

4>n ■= log hn o J, and it converges to 4> := log ho J. As we saw in the previous section, it 
follows that {(C0„)±} converge uniformly on T to (C0)*, respectively. Thus, the claim 
follows. 

Finally, we compute the Szego function of h{a,x; •), a G (0, 1], a; G -F, defined in 



(2.161. It is a simple observation that 



2(z — x) 

= 1 and S'n(i,i:;.)(z) = — z e D, 

by uniqueness of the Szego function. Let us also recall that S^i ^-..-^ is a non-vanishing 
holomorphic function in D, including at infinity, with zero winding number on any curve 
there, and = 1- Thus, for any a G (0, 1], an a-th power of 6*^1(1.2;:.) exists. In 

other words, the function 

Sa ^h(l,x;-)^ ^a{oo) — 1, 

is a single-valued holomorphic function in D with continuous traces on each side of F. 

Denote by the branch of the argument of (p which is continuous in C \ {(— cxd, y] U 
Fgu x^F} and vanishes on positive reals greater than 1, where Faux was defined in Section 
2.3 Then 9^ has continuous extensions to F that satisfy 0+ = —9~ on F, 9^{1) = 0, 
and 9^{—l) — tt. It can be easily checked that AigS^ — a{9x — 9^), where AvgSa is the 
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continuous branch of the argument of Sa that vanishes at infinity and 9x was defined in 
Section [2^ Thus, 

± \2{t-xT \ exp{ia(0,(i)-0±(t))}, t e 

^ |(p±(t)|" \ exp{*a(0±(i)-0±(i))}, teF_i,,\{x}, 



and therefore S'^S~ = h{a,x; ■) on F by (2.6 1 and since 9^ = 6*+ — 27r there. Altogether, 
we have proved that Sn(a,x;-) = Sa- 

Now, define Cq, := / on i^. Then, for x ^ ±1, we have 

r (f) - |^-/.^|2a / exp{-ia(0+(i) 

- ^'^1 \ exp{~^ai-etit) + 0-{t)+e+{t)-9-{t))} 

^ ^' \ exp{-2ia6'+(i) + 2ia7r}, t e \ {x} ' 

for a; = 1, we have that 

Ca{t) = \(p-{t)f'^exp{~2ia0+{t) + 2iaTT} , t e F; 

and finally, for x — —1, we have that 

Cait) = |^-(i)p"exp{-2ia0+(t)} , teF 

Therefore, if a; 7^ ±1 then and are piece wise continuous functions on F, having 
jump-type singularities at x of magnitudes 2 sin(a7r) /(^+(a;)|" and 2 sin(a7r)|(y9+(x) 

respectively, and Cq(±1) = 1; if a; is ±1 then Cq is a continuous function on F, Cq(^1) = 1, 



and Cq(±1) = exp{±2iQ;7r}. Thus, (2.17) does indeed hold. 

Now, let Fq C F be a finite set of distinct points. Given a,j. < 1 for every x G Fg, 
let h be associated to Fq via ( |2.18[ ). Then 

Qa;/2,a:;-) 

a:ei="o\{±l} xGFon{±l} 

by the multiplicativity property of Szego functions. Moreover, = / Sf^ and = 
Shl^h ^'"^ products of Ca^ and 1/cq^, respectively, over a; G Fq, and therefore they 
are piecewise continuous function on F with jump-type singularities at cc G Fq \ {±1}. 

3.4. Smirnov Classes 

Let r be an analytic Jordan curve. Denote by and the bounded and unbounded 
components of C \ F. Suppose also that T is oriented counter-clockwise, i.e. £>+ lies on 
the left when Y is traversed in positive direction. Assume also that G . Set U", 
p € [IjCxd), to be the space of p-summable functions on T with respect to arclength and 
denote the corresponding norms by || • \\p. For each cj) E the Cauchy integral operator 



( 3.2 ) defines a sectionally holomorphic function in UD . To be more specific we shall 
denote 

C±0=(C0)p±. 

Then each function C^4> is holomorphic in the corresponding domain and has non- 
tangential boundary values a.e. on F, which, as before, we denote by {C(f>)^. For each 
p € [1, 00) we define the Smirnov classes F^ and F^ as 

El := {(/.eFP; (f, = O} , 

EP_ {0eFf : C+0 = (C+0)(O)} . 
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Clearly, Ej. are closed subspaces of L^. We note also that the Cauchy formula holds, i.e. 
E^ then 

{C(l>)+ = (l> and (Ccf))' = -(p+{C+4>){0) a.e. on T 



by Sokhotski-Plemelj formula (3.3 1. In other words, (j) admits an analytic extension to 
whose L^-means turn out to be uniformly bounded on the level curves of any conformal 
map from D onto D^. Conversely, if /i is a holomorphic function in whose L^- 
means are uniformly bounded on some sequence of rectifiable Jordan curves in whose 
bounded or unbounded components of the complement exhaust or D~ , respectively, 
then the trace of h belongs to E^. (see |TB] and [TT] for details) . 

It will be of further convenience to define the projection operators from to the 
Smirnov classes E^. To do that, we first define the Cauchy singular operator on F, say 
S, by the rule 



{sm f ^dr, t e r, 

m Jy T — t 



where the integral is understood in the sense of the principal value. We also set I to 
stand for the identity operator on and we put E^, p G [1, oo), to denote the subspace 
of E^ consisting of functions vanishing at infinity, i.e. 

EP_ {c^^E'L: (C+(/))(0) = 0} = {0 G : C+0 = O} . 

In the considerations below we always assume that p G (1, oo). Then 

LP = El® El. (3.11) 

We define complementary projections on L^ in the following manner: 

V: LP^El, 0^(l/2)(I+5)0, 

Q: LP^El, <j>>^{l/2){I-S)cf>. ^''■'^> 

These operators are bounded and can be used to rewrite the Sokhotski-Plemelj formula 
as 

rcl)={C(l))+ and Qcl} = -{C(j))-, (j) & L^ . (3.13) 

Finally, we characterize the weak convergence in El. According to the Riesz rep- 
resentation theorem, any linear functional £ G {L^)* has a unique representation of the 
form 

£cl)= f cbf\dT\, cbeEP^, feL\ -+^=1. (3.14) 
Jv P Q 

It follows now from the Hanh-Banach theorem that 

{El_r - Ly{EP+)^, (3.15) 

where the symbol "=" stands for "isometrically isomorphic" and (El)-^ is the annihilator 



of El under the pairing (3.141. It can be easily checked that 
Since F is analytic, multiplication by dr/ldrj is an isometric isomorphism of L"^ into itself. 



: l^eEl]. (3.16) 



Thus, by (3.11 1 and (3.121 every f £ L"^ can be uniquely decomposed as 



18 



L. Baratchart and M. Yattselev 



Combining (3.151, (3.16), (3.17), and (3.11), we see that for any linear functional t £ 



{E^)* there uniquely exists ip £ E'L such that 

= y (j)iPdT, (j)£El. 

Since (• — z)^^ e E'i_ for every z G D'^ , it follows from the Cauchy formula that a bounded 
sequence {4>n\ C E^ weakly converges to zero if and only if {C+0„} converges to zero 
locally uniformly in Z?^". 

3.5. Toeplitz and Hankel Operators 

In the notation of the previous section, let c <E L°° , i.e. c be a bounded function on T 
with the norm ||c||oo, and let p € (1,cxd). We define the Toeplitz and Hankel operators 
with symbol c by the rule 

%: El^El, <^^P(c0), 
He-. El^El, 4,>^Q{c^), 

respectively. Then 7^ and He are bounded operators and He is compact whenever c G 
C(r) (see, for instance, HH Prop. 6.21]). 

Let {(j)n} C E^ be a bounded sequence that weakly converges to zero and {c„} be 
a relatively compact family in L°° . Then for any limit point of {c„}, say c, we have 

||Tc(/)„ ~Tc^(l)n\\p = ||7'((C- Cn)(f>n)\\p < ||7'||||c- C„||oo ||0n||p, 

where ||7^|j is the operator norm of V. Therefore, since V and 7^ are bounded operators, 
it can be readily verified that {7^^(/)„} weakly converges to zero. Moreover, using com- 
pactness of Hankel operators with continuous symbols, a similar argument shows that 
{Hc„(j)n} strongly converges to zero in the || • ||p-norm if {c„} is a relatively compact 
family in C(T). 

Let now c be a continuous function on F except for a finite set of points, say Tq. 
Assume further that c possesses finite one-sided limits, say c(C^) and c(C~), at each 
C G Tq. Then there exists c' € C(T) such that 

l|c-c'||oo = Jmaxj,(C), MO ■■= |c(C+) -c(r)|. 

Thus, if {4>n} C E^, \\<f)n\\ — 1, is weakly convergent to zero, it holds that 
limsup ||7Yc0n||p < limsup(||Hc-c'</'n||p + ||Wc'0«||p) 

n — >-oo 71 — >oo 

< lim sup II Q((c - c')0n) lip < ^IIQII maxjc(C) 
since He' is a compact operator. 



3.6. On Conjugate Functions 

Let r be an analytic Jordan curve and and stand, as before, for the bounded and 
unbounded components of the complement of F. Denote by / a conformal map of the 
unit disk, D, onto D~^. Further, let be a sequence of holomorphic functions in D+ 

that continuously extend to F. Assume also that {|(/>n|} is a relatively compact family in 
DC*{T). Define ■0n ■— 't'nof- Then {■0n} is a sequence of bounded holomorphic functions 
in D with continuous boundary values on the unit circle, T. Moreover, since /' extends 
continuously on T by HSl Thm. 3.3.5], it holds that {iV'rtl} is a relatively compact family 
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in Z?C* (T) (as before, we use the same symbol ui since only the existence of the majorant 
is important). 

By the canonical factorization of Hardy functions in D (see [111 Sec. 2.4]), we get 

ipn = 4'ni^n7 where IV'^I = 1 on T and 

C(^):=exp|^^log|^„(T)|^|, ZED. 

In other words, "0^ and t/'J^ are the inner and outer factors of ipn, respectively. Since 
\ip°\ = \i/jn\ on T, {IV'nl} is a relatively compact family in DC^{T). Moreover, since \ipn\ 
is uniformly bounded below on T, we have that 

^iog\i,-\{S) < cj|^o|((5)/min|V'°(r)| < kiUj{S) 

r t Jl 

for all S small enough, where ki is a positive constant independent of n. Thus, without loss 
of generality we may assume that {log IV'nll ^ DCuj{T). Denote by argi/;" the harmonic 
conjugate of log normalized to be zero at zero. Then by [19, Thm. III. 1.3] we have 
that 

where k2 is an absolute constant and the moduli of continuity of arg are computed on 
T. This, in particular, implies that {arg?/;°} is a uniformly equicontinuous family of uni- 
formly bounded functions on T. Hence, by Ascoli-Arzela theorem, {argtp"} is relatively 
compact in C(T) and therefore so is {"iAnl- 

Summarizing the above we arrive at the following. If T is an analytic Jordan curve 
and {(j)n} is a sequence of holomorphic functions in with continuous traces on T such 
that {|0n|} forms a relatively compact family in DC*(T), then each 0„ can be written 
as 

K^^nK, \M^m on r, 

and the sequence {0^} forms a relatively compact family in C(T). 



4. Proofs 



In the first part of this section we prove Theorems |2] [3] and |4] A number of places in 
these proofs contain references to the material of the previous section. The second part 
is devoted to the proof of Theorem [T] and is independent of the rest of the paper. 

4.1. Proofs of Theorems [2||3} and |4] 

Proof of Theorem The starting point for the proof of this theorem is the method of 
singular integral equation proposed by J. Nuttall in |38| . 

By the very definition of (3.6), (2.6l, and (2.121, it follows that 



(i?„w)+ + (i?„w) = 2qnWn = 2G^^q„ 



The last equation can be rewritten as 



= 2G,. 



on F. 



QnSu 



(4.1) 



(4.2) 
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on F. Observe that 



(4.3) 



is a holoniorphic function in D and G C{F), a„ (±1) — a^(±l). Indeed, Sw„ is a 
non- vanishing and holoniorphic in D, 

w(z)^"(z) = 0(z"+i) as z->oo 

by the very definition of these functions, and 

i?„(z) = 0(l/z"+i) as z^cx) 



due to orthogonahty relations (2.2 1. The continuity of the boundary values and equal- 
ity of at ±1 follow from (3.6) and (3.9). Let also 



2Gu 



If" 



(4.4) 



It is again immediate that 6„ is a holomorphic function in D, G C!{F), and &^(±1) = 
6~(±1). Then using the multiplicativity property of Szego functions, we can rewrite (4.2 1 
as 



+ {anCnY 



where 



= 5J on F, 



F. 



(4.5) 



(4.6) 



Clearly, e C{F), c^c^^ = 1, and c+(±l) — c„ (±1) = 1. Moreover, {c^} are relatively 
compact families in C{F) (see Section 3.3 the discussion right after (3.10)). Further, by 
( |2.12| , p!6| ), and ( |3.10| it holds that 

(^"5^J± 1 i;„ 



7'J on F. 



Thus, we obtain from (4.5) that 

6? 



+ {anTnCn)^ On F. 



(4.7) 
(4.8) 

J^^{F) is an analytic 
curve (see Section 3.1 1. Denote by £)+ and the bounded and unbounded components 
of the complement of F, respectively (we orient F counter-clockwise), and define 

f-^f°J\D+, f^a.a,bn,rn. 

Then a„, b„, f„ are holomorphic functions in Z?+ with continuous traces 



Let J be the usual Joukovski transformation. Then F 



(4.9) 



Continuity of the boundary values immediately implies that these functions are bounded 
in D+ and therefore 



feEl and /(l/-)ei?^, f 



(4.10) 



for any p e [1, oo) by the characterization of Smirnov classes (Section 3.4 1. Moreover, the 
sequence {r„} converges to zero locally uniformly in _D+ and therefore weakly goes to 
zero in E^ (see Section |3.4| by the very definition of f„ and the fact that F is symmetric 
with respect to S. It follows from (14. 81) that 



6„(r) = a„(l/r) (a„f„c„)(r), r G F, 



(4.11) 
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where c„|r± '■— o J. Notice that c„ are, in fact, continuous functions on T and {c„} is 
a relatively compact family in C{T). Fix an arbitrary p G (1, oo) and denote 

on £)+ and rf„,p(-) := a„,p(0) - a„,p(l/-) on D", (4.12) 



■*"n \ \p 



where || • ||p is the L^-norm on F. Clearly, d^^p G E"^. Applying subsequently the operators 



V and Q (see ( [3l2| )) to both sides of ( [4TT| , we get in the view of ( [4lO| that 



(4.13) 



Equivalently, we can rewrite (4.131 as 



(4.14) 



^n/||^n||p — ^n,p(0) ^^(^n,p^n) 

where 7^ and 7Y„ are the Toeplitz and Hankel operators with symbol c„, respectively 
(see Section 3.51. By the normalization of {a„ p} and the properties of {r„}, the norms 
IjO'n.p'^nllp are uniformly bounded with n and the sequence {a„^pf„} weakly goes to zero. 
Thus, 

IMn,p||p ~* as n — > CX3 

by the compactness properties of Hankel operators with continuous symbols (see discus- 
sion in Section 3.5 1. In particular, since 

l!<p(i/-)llp < 



mm^gr \t 



2/p 



as n 



(note that T cannot pass through the origin since otherwise it would be unbounded as r 
and 1/r belong to T simultaneously), we have from (4.121 that 



|a„(0)|/||a„||p 1 
||1 - a„/a„(0)||p 



oo. 



(4.15) 



Using now properties of Toeplitz operators with L°° symbols (see Section 3.5 1, we derive 
from (4.141, (4.151, and the Cauchy formula that 



o(l) 
-o(l) 



(4.16) 



bn/an{0) = 
an/an(0) = 

locally uniformly in or equivalently 

a„(0) ^" 

1 RnWip"^ 

. a„(0) Sio„ 

locally uniformly in D. The first asymptotic formula in (4.16) immediately implies that 
Qn has exact degree n for all n large enough and therefore can be normalized to be 
monic. Under such a normalization we deduce again from the first formula in (4.16) by 
considering the point at infinity and using the normalization Su)„(oo) — 1 that 

a„(0) = [1 + o(l)]2i-"G^„ = [1 + o(l)]2"7„ (4.17) 



o(l) 
o(l) 



and (2.13) follows. 
Now, let 



RnW 
^nSn 



2"7„ 



(4.18) 
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Then dn '■— dn o J|c+ is such that 

IMnllp < ||1 - an/a„(0)||p + o(l)||a„||p/|a„(0)| ^ as n -> oo 



(4.19) 



for any p e (1, cx)) by (4.17) and (4.15 1, and therefore also for p = 1 by Holder's inequality. 



So, we have for any p G [1, oo) that 

\dn{rW\dT\ 



{\d-{t)n^^{t)Y\ + \d+{t)n^-{t)Y\)\dt\ 



{\d-{m^^{t)\ + \d+{t)n^-{t)\) 



dt 



w+{t) 



(4.20) 



as n ^ oo. Since \(p^\ > minr |t| and |w"'"(t)| — — t'^\, the second part of (2.141 
follows from (4.181 and (4.201. Moreover, since we can rewrite (4.11 as 



g„ InS+S- = {R„w)+ + (Rnw)- = 7„ [((1 + d„)Sn)+ + ((1 + d„)5„)- 



by (4.181, we get the first part of (2.141 



So, it remains only to prove (2.15 1. By Lavrentiev's theorem, which is the analog of 



Weierstrass theorem for complex-valued functions on an arc, it is enough to show that 



dt 



w+{t) 



F W+{t) 



, J = 0,1,. 



(4.21) 



By (4.11, we get 



47„W„ 



{{Rnw)+y + 2{R„w)+{Rnw)- + {{Rnw)-y 



In' 



l{{RnW) + f , {RnW)+{RnW)- , I {{RnW)-f 



2 Sn Sn Sn Sn ^ Sn Sn 

Recalling the definitions of a„, r„, c„, and d„, we get in view of ( |4.7| that 

In^qlwn = 2"2"7;^^ [(a^,r„c„) + /2 + a+a,7 + (a^r„c„)"/2] 
= 1 + [d- + d+ + d-d+] 

+ [((1 + d„)2r„c„)+ + ((1 + d„)2r„c„)-] /2. 



Thus, 



where 



and 



7„w+(i) 



w+{t) 



h,n — T^ih,n, 



Pdt 



h,n [d-{t) + dtAt) + d-{t)dt{t)] ^^^^^ 



l2,n — 



2m 



Pdt 



[(1 + dl{t))\rU-^){t) + (1 + d-{t)f{r-c-){t)\ — . 



Equations (4.211 shall follow upon showing that /i^„ and /2,n converge to zero as n 
grows large. Clearly, this holds for /i_„ by (4.20) since are bounded on F. Further, 
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substituting t ~ J(t), t e F, we get as in (3.5) that 



^(l + d±(t))^(r„c„)±(i)^ 



(l + d„(T))2(f„c„)(r)J^"(T) 



dr 



FT 



where we integrate on T from —1 to 1 and on r+ from 1 to —1. Thus, by the very 
definition of the operators C and V (see (3.2 1 and ( |3.12 1), it holds that 



;C+((l+d„)2f„C„J^)j (0) - (C+P((l + d„)'^nC„J^)) (0) 

= (C+7^-„((l + d„)2f„)) (0), 



where 7^-.„ 



is the Toeplitz operator with symbol c„J^. Since ||1 + fi„||p = 1 + o(l) for 

are uniformly bounded above on T while {r„} 



any p G (l,oo) by (4.191 and since |f 
converges to zero locally uniformly in Z)+ by the assumptions of the theorem, the sequence 
{(1 + dn)'^fn} weakly goes to zero in each L^, p e (1,cxd). Furthermore, {c„J-'} is a 
relatively compact family in C(r) for every j — 0,1,.... Therefore TJ^„((1 + (in)^f„) 
converges weakly to zero in any L^, p g (1, oo), and subsequently /2.n tend to zero as n 
grows large (see Section |3.5[) . This finishes the proof of the theorem. □ 



Proof of Theorem^ The following is a modification of the proof of Theorem [2] Thus, 
we shall keep all the notation we used in that proof. 
Set C|r± 



o J, where 



to (4.11 ), we get that 



were defined after (2.18). Repeating the steps leading 



6„(t) = a„(l/r) + {anfnCnc){T), r e r. 



(4.22) 



As before, 6„ and f„ are holomorphic functions in with continuous traces on F, {r„} 
converges to zero locally uniformly in and the boundary values of |f„| are uniformly 
bounded above on F. Further, {c„} is a relatively compact family in C(F) and a„ are 
holomorphic functions in _D+. However, the traces of a„ are continuous only on F \ Fq, 
Fo = J^^(Fo). Moreover, at each C G Fq, each a„ has an algebraic singularity of order 
—ax > —1/2, X — J{C). Thus, by the characterization of Smirnov classes (Section 3.4), 
we have that a„ G E'^. Finally, the function c is piecewise continuous on F, i.e. c is 
continuous on F \ Fq with jump-type discontinuities at each C e Fg. 



Given (4.22), we obtain as in (4.13) that 



bn/Wdnh = a™, 2(0) + 'P(a,i,2^nC„c) 
dn,2 = Q{dn,2fnCnC) 



(4.23) 



By assumptions on S', it holds that {|f„|} is a relatively compact family in DC^{T). 
Thus, as explained in Section 3.6 we have that f„ — fl^fn, where {f°} is a normal family 



we have that fn — f 

in £)"*" with continuous boundary values on F that form a relatively compact family in 
C(F). Moreover, |f„| = |f°| on F and since neither of the limit points of {f°} can vanish 
in by the outer character of f°, it follows from the assumptions on r„ that {fjj} 
converges to zero locally uniformly in and |f*J = 1 on F. 



Now, we can write (4.23) as 



^n/llOnlb = an,2(;0) - T„(a„^2?^n) 
dn,2 = 'Hn{dn,2rn) 



(4.24) 
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where Tn is the Toephtz operator with symbol c„c and 7i„ is the Hankel operator with 
symbol c„f°c. Then we get from the discussion at the end of Section [33] that 

limsup ||(i„,2||2 < ^IIQII maxlimsup |(c„f°c)(C'*") - (c„f°c)(C")| 

since {c„f°} is a relatively compact family in C{T), ||an,2'^^||2 = 1, and {fin, 2?^^} weakly 
converges to zero in i?^. Then it immediately follows from the very definition of c„, c, 
f„, and ( |2.19| ) that 

limsup ||d„,2(l/-)ll2 < 1 

n — >oo 

and since (i„.2(l/-) = (an(0) — an('))/ll^n||2i we derive that 

limsup ||a„||2/|a„(0)| < oo. (4.25) 

n — ^oo 

As Tn and 7i„ are bounded operators acting on weakly convergent sequences, we get 
from (|4.24|) and (|4.25| that 



6„/a„(0) - l + o(l) 
a«/a«(0) = l + o(l) 



(4.26) 



locally uniformly in D+. As in the proof of Theorem [2] ( |4.26 1 implies (4.16) and (4.171 
and therefore (2.13) follows under the present assumptions. 



Define dn as in (4.18). Then it still holds that 

dn = dn o J|£,+ = [1 + 0(1)] :r-^ - 1 
and hence limsup„_^g^ ll'^nlb < 00 by (4.25). Thus, as in (4.20), we get 



limsup / {\d-{t)\'\^+{t)\ + \d+{t)\^\^-{t)\) 



w+{t) 



The rest of the conclusions of the theorem now follow as in the proof of Theorem [2] from 
the above estimate. □ 

Proof of Theorem^ The existence of schemes that make F symmetric with respect to 
them was shown in Theorem [T] It is also well-known and follows easily from the defining 
properties of Pade approximants, that the denominators of n„ satisfy non-Hermitian 
orthogonality relation of the form 



t=qn{t)^^0, J = 0, 
Vn(tj 



and the error of approximation is given by 



(^-n„)(z) 



Vn{z) f ql{t) dii{t) 



Qniz) J Vn{t) Z-t 



z £ C \ supp(/^). 



Given (2.21), we see that the asymptotic behavior of g„ is governed by (2.13). Thus, 
using the orthogonality of qn, it is not difficult to see that we can write 



in _D = C \ supp(/i). 
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where i?„ is the function of the second kind associated to g„ via (2.3 1. Hence, equations 



(2.131 imply that 



(1 + 0(1)) 



^nVnSl 2G^+o(l) Si 



w 



2Ga + o(1) ^2 ^ 



locally uniformly in D, where we used (3.10) 



□ 



4.2. Proof of Theorem [T] 

Before we proceed with the main part of the proof, we claim that our assumptions on F 
yield the existence of the partial derivatives dU^~^ / dri^ almost everywhere on F, where 
A is the weighted equilibrium distribution on F in the field —V^ and v \s a, compactly 
supported probability Borel measure in D = <C \ F . 

Let '0 be the conformal map of D into the unit disk, D, such that ij^ipo) — and 
^''(oo) > 0. Put 

^{z) - V'(e) 



5(e;z) := 



e,z e D. 



(4.27) 



l-7A(e)^(z)' 

Then 6(e; •) is a holomorphic function in D with a simple zero at e and unimodular 
boundary values on F. Hence, the Green Junction for D with pole at e € D j49| Sec. II. 4] 
is given by — log |6(e; •)! and the Green potential of a Borel measure v, supp(j^) C D, not 
necessarily compact, is given by 

U^iz) :=- j \og\h{t-z)\dv{t), zeD. (4.28) 

Furthermore, since D is simply connected and therefore regular with respect to the 
Dirichlet problem, extends continuously to F and is identically zero there. 

Assume now that v has compact support. Then it is an easy consequence of the 
characterization of the weighted equilibrium distribution [49L Thm. 1.3.3] and the repre- 
sentation of Green potentials via logarithmic potentials [49l Thm. II. 5.1] that A as above 
is, in fact, the balayage [49j Sec. II. 4] of ly onto F, supp(A) = F, and 



{U^ - U''){z) ^ -U^aiz) - log\b{^;t)\d,y{t), zE 



(4.29) 



Thus, the existence of dU^ ^ /dri^ specializes to the existence of dlfj^/dn^. 

Let now J be the Joukovski transformation and F J^^{F). It was explained in 



Section 3.1 that F is a Jordan curve. Denote by Z)+ and D the bounded and unbounded 



components of the complement of F, respectively. Then 



(I) -.^ o J, : D^ 
is a conformal map. Set 

u{z) :--(C/^o J)(z) = 



0(0) = 0, and 0'(O) > 0, 



log 



(z) - m 



0(z)0(t) 



(4.30) 



where di'~^ :— d{vo J\y)+). As J is conformal on F \ {±1}, the existence of du/dn a.e. on 
F, where d/dn is a partial derivative with respect to the inner normal on F, will imply 
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the existence of dUj^/dn^ a.e. on F. As 

Ou f dvu 

-— (r) lim (n,-, Vu(r + Sur)) = 2 lim Re nr-^—ir + Jn,-) 

where (•, •) is the usual scalar product in M^, Ut is the inner normal vector at r, rir is the 
unimodular complex number corresponding to fir, and (d/dz) :— {{d/dx) — i{d/dy))/2, 
it is needed to show that the function du/dz, holomorphic in \ supp(:^+), has non- 
tangential limits a.e. on T. It can be readily verified that 

9u mf^^^mi ,,n-. (4.31) 

dz' ' 2 y (l-0(r)0(z))(,^(z)-</.(r)) 

for the definition of 



3.4 



Moreover, by construction, we have that (j>' G E]_ (see Section 
the Smirnov class E^). Recall (see Section 3.1 1 also that the assumed condition on the 
behavior of F near the endpoints implies rectifiability of F. Hence, any function in E\_ 
has non-tangential limits a.e. on F by [TSJ Thm. 10.3]. As (j) extends continuously on F 
by Caratheodory's theorem [45l Sec. 2.1], the desired claim on existence of the boundary 
values of du/dz and respectively normal derivatives dU^"" / dn^ follows. 
(a)^(b): Let F be symmetric with respect to a triangular scheme of points S" = {En}- 
Denote by i/„ the counting measures of points in En and let be a weak* limit point of 
{t^n}. In other words, there exists a subsequence Ni C N such that ^ z/, Ni 9 n — > oo. 
If ly has bounded support, we shall show that F is symmetric in the field —U" . By (2.9 1 
and (4.291, this is equivalent to proving that 

If v has unbounded support, we claim that there exists a compactly supported Borel 
measure v* such that = in some open neighborhood of F. Thus, the equality in 
(4.321 will be sufficient to show that F is symmetric in the field —W^ . To construct i^*, 
set k{z) = l/{z — a), a ^ supp(z^) U F. Clearly, = ^^^(13) ° where di^i^ = d{v o k~^) 
has compact support. Let i be a Jordan curve in k{D) encompassing supp(i'k). Then the 
balayage of on L, say i^* , is such that U^^j^-^ = ^k.{d) outside of L [IS] Thm. II.4.1]. 
Thus, the claim follows for d{iy*) := d{i'* o k). 

In order to prove that the equality in (4.321 indeed holds, we derive different from 



(4.281 integral representation for the Green potential of the above measure v. It follows 
from ( |4.28[ ) that 

C/^"(2) = -^log|6„(z)l, bn{z) := Y[ b{e;z), z e D, neNi. 

Further, since F is symmetric with respect to S', {log lr„/6„]} is a sequence of harmonic 
functions in D whose boundary values uniformly bounded above and away from zero. 
Hence, we get from the maximum principle for harmonic functions that 

-^(loglr„(z)l-logl5„(z)l) ^0 uniformly in C. (4.33) 
Zn 

Moreover, removing o{n) factors simultaneously from r„ and 6„ will not alter this con- 
clusion. Indeed, new difference will remain harmonic and its boundary values on F will 
still uniformly go to zero since |&^(i; ■)! = 1 ^^d |r^(t; •)] < const, for all t G supp(^), 
where const, depends only on supp((f). On the other hand, since f„ weakly converge 
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to u, it is possible to write each Vn as a sum of two measures, say z/„^i and 1^,1,2, such 
that \i'n,2 \ — o{n) and the supports of Vn,i are asymptotically contained in any open set 
around supp(i/). Let r„ 1 and 6„ 1 correspond to Vn.i as r„ and 6„ correspond to v^. Then 
we obtain from the weak* convergence of Un.i to v and the remark after (4.331 that 



lim Uj^ 



lim 

Ni 9n — »-oo 



1 

2n 



log |r„ 



locally uniformly in _D \ supp(i^). So, writing (l/2n) log |r„_i | as an integral of log |r(-; z)| 
against Vn,!, we derive from the weak* convergence of measures that 



\og\r{t; z)\dv{t), z e £> \ supp(z^). 



(4.34) 



We proceed by reformulating (4.32 1 in terms of the boundary values of the complex- 
valued function H dV^/dz. Since is harmonic in £> \ supp(i^), H is holomorphic 
there on account of Cauchy-Riemann equations [T, Sec. 4.6.1]. Moreover, it follows from 
\AM\ that 



H{z) 



^'{z) 



-dv{T), 



(4.35) 



(1 - Lp{T)ip{z)){ip{z) - Ip{t)) 

z g D\supp(z/). Observe further that H extends analytically across each side of F\{±1}. 
(Note that ip' = Lp/w and therefore it extends analytically across each side of \ {±1} 
as and w obviously do.) 

Let now t e F^, where Fr C F \ {±1} is the set of points at which F possesses 

tangents. Denote by ft and the unit tangent vector and the one-sided unit normal 
vectors to F at t. Further, put and nf^ to be the corresponding unimodular complex 
numbers, — irt. Then 

(VC/^(z), n±) = 2Re (n± H{z)) , z e D. 

As H extends holomorphically across F \ {±1}, the above equality also holds at z — t. 
In other words, we have that 



-^it) = 2Re {nfH^{t)) = 2n±iJ±(t), t e F,. 



(4.36) 



The last equality in (4.361 is valid because = on F and therefore 



= ^(t) = T2Im {nfH^{t)) , t e F^ 



As nl 



If. , (4.32) will follow from (4.361 if we show that 

H+ = -H- 



on Fr 



(4.37) 



The latter is well-understood in the theory of symmetric contours [24, pg. 335] and can 
be seen as follows. 

Observe that (<y9±)' = ((/?')± on F \ {±1} since ip extends holomorphically across 
each side of F\ {±1}. Therefore, we have by (2.6) that 



w (p ) w+ 



on F\{±1}. 
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Then by (4.35 1 we have for the unrestricted boundary values of H on each side of F that 



2 

2(^+)2(t) 



(l-<^(r)^+(i))(^+(t) 
1 - ^^r) 



V{r)) 



dv{T) 



(^-(t)-^(T))(l-^(r)^-(i)) 
which finishes this part of the proof. 

(b)=J>(c): Let F be symmetric in the field —V^, where v is a, positive Borel measure 
compactly supported in D. We show that there exists a univalent function p holomorphic 
in some neighborhood of [—1, 1] such that F = p([— 1, 1]). 

Let r, and be as before and denote by the domain in C such that 
J{n) = C\supp(i/). Set 

Then u is identically zero on F, harmonic and negative in , and harmonic and positive 
in D~ . Moreover, u^d+ has integral representation (4.301 and, as explained after (4.311, 
it has well-defined non-tangential boundary values, say {du/dz)~^, that are Integra 
r. Since u{l/z) = —u{z) and = {1/z : z E D^}, the trace {du/dz 
belongs to L^. 

In another connection, it is easy to verify that 

{du/dz){z) ^TH{J{z))J'{z), zenoD"^, H:^dU^/dz, (4.38) 



Die on 
also exists and 



is a sectionally holomorphic function. Let ip be given by (2.5 1. Then F is the boundary of 
ip{D) and, as before, we set F* = (p^{F). Since (2.9) is equivalent to (4.371 as explained 
in the previous part of the proof, we obtain for a.e. r € F~ that 

{du/dz)+ (t) = -H+{J{t))J'{t) = H'{J{t))J'{t) = (du/dz)- (r). 

Analogously, we can show that (du/dz)^ coinsides with (du/dz)^ a.e. on F+ and there- 
fore a.e. on F. It follows now from an application of the Cauchy formula that du/dz is 
analytic across F, i.e. du/dz is a holomorphic function in the whole fi. 

Recall that u is identically zero on F and does not vanish anywhere else in C. Hence, 
the level lines F^ := {z : u{z) = e} are single Jordan curves for all e sufficiently close 
to zero. Let L be such a level line in and fii be the annular domain bounded by 
L and L^^ :— {z : 1/z e L}. Observed that by Sard's theorem on regular values of 
smooth functions we always can assume that L is smooth. Since u is constant on L, the 
tangential derivative of u there is zero and we get as in (4.36) that 

du 



du 

^^-(r) - (Vu(r),n,) = 2rir^(T), 



T e L, 



(4.39) 



where d/dn is the derivative in the direction of the inner normal with respect to fl^, 
and rir is the corresponding unimodular complex number. As —u is the Green potential 
of a probability measure v'^ (see (4.30 1), L is smooth, and dr = irirds on L, it follows 
from ( |4.39[ ) and Gauss' theorem [49, Thm. II.l.l] that 

du 



2m = 27rw+{D+) = i 



dn 



ds^2 



du 
dz 



{T)dT. 



(4.40) 



Furthermore, established for L, the chain of equalities above holds, in fact, for any Jordan 
curve contained in fi^ and homologous to L since du/dz is analytic on fi^, where we 
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always take the outer normal with respect to the inner domain of that Jordan curve. 
Thus, a function 

du 



■^{z) :=exp<i 2 j g^Mdr 



(4.41) 



is well-defined and holomorphic in Cl^, where the integral is taken over any path joining 
1 and z and lying entirely in il^. 

Now, observe that the maximum principle for harmonic functions yields that du/dn 
is non-negative on each contained in fi^. Set 7e(T) to be the subarc of obtained by 
traversing it into the counter-clockwise direction from to t, where {x^} := D (0, oo). 
Then ^^^{du/ dn)ds is a positive strictly increasing function of r with the range [0, 27r] 
as r winds once on in the positive direction starting from 1 by ( |4.40[ ) and the remark 
right after. Hence, choosing the initial path of integration to be 7(t) [^^Xf] U 7£(t), 
we derive as in (4.40) that 



logl^'(T) 



Re 



= 2Re 




du 




du 
dn 



ds 



du , 



(C)dC =:logp„ 



and 



Arg('i'(T)) = Im 2 



hi-r) 



du 
dz 



(C)dC = 2Im 



du 
dz 



7.(r) 



dn 



ds, 



r € Fg, where Arg(^') is the principal value of the argument of Hence, ^> is univalent 
on each F^ C f^L and *(F,) = Tp, (in particular, ^'(F) T and ^'(±1) = ±1). As each 
level line F^ lies either entirely inside of J^l or entirely outside, ^' is a univalent function 
in the whole domain fi^. It also holds that 



^{l/z)^l/^{z), zeilL, 



(4.42) 



as follows from a change of variables in (4.41 1 and since 

{du/dz){l/z) = ±H{J{l/z))J'{l/z) ^ Tz^H{J{z))J'{z) = z^{du/dz){z), 



zeflLnD"^, by (4.381 



Set / := JovJ/ot^. Since ip is univalent and holomorphic in J(rii)\F, ip{J(flL)\F) — 
UlD , and is injective and holomorphic in il^, / is also univalent and holomorphic 
in J{nL) \ F. As ip^{F) = F± and «'(F±) = T±, it holds that /+ and /- both map F 
onto [—1,1]. It is also true that 

nt) = j(vi/(^-(t))) = j(vi/(i/^+(<))) = j(i/vi/(^+(t))) 

= J(vi/((^+(<))) = /+(<), teF, 



by (2.61 and (4.42). Hence, / is a holomorphic and univalent function in some J{flL) 
that maps F onto [—1,1] and therefore the desired analytic parametrization of F is given 
by P = /"^ 

(c)=>(a): Let F be an analytic Jordan arc and p be its holomorphic univalent parametriza- 
tion. Denote, as usual, by F the preimage of F under the Joukovski transformation J. It 
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was shown in Section 3.1 that there exists a function ^, holomorphic in some neighbor- 
hood of r, say n, such that ^'(r) = T and 

^-(1/2) = l/^'(z), z,l/z€n. (4.43) 



In fact, $ = ^ was constructed in Section 3.1 Let p € (0, 1) be such that 

Cp(0),Ci/p(0) C^'(r!), a,izo):={zeC: \z-zo\^x}, x > 0. 

Denote Tp $(Cp(0)) and T^/p := $(Ci/p(0)). It immediately follows from ( |4.43| that 

ri/p-r;i = {reC: l/reTp}. (4.44) 

Denote by fip the annular domain bounded by Tp and r^/p and define 

u{z) := log |^'(2;)|, z e f^. 

Then u is a harmonic function in some neighborhood of flp such that u = on F, u = log p 
on Tp, and u = — log p on Ti/p. Furthermore, it follows right away from the definition of 
u that 

du, . l^'{z) 



(4.45) 



Let now rir stand for the unimodular complex number corresponding to the inner normal 
of rip at r g dQ.p. Since inner normals of 9^'(rip) are represented by ^/|^|, ^ G C'p(O), 
and — C/ICI, £, € Ci/p(0), we have by conformality at ^ that 



± 



|M/(r)| |M/'(t) 



r G Fpii. 



(4.46) 



As u is harmonic in a neighborhood of fip, we deduce from (4.45) and (4.461 that 
^""(r) = (Vu(r),nO = 2Re (^{T)n^ = ±p^i|vl/'(r)|, r G Fp±i, 



9n ' \9z 

where n-r is the unit vector corresponding to the complex number n-r. Observe that 

^-'(1/2) _ ^'{z) 

by ( |4.43| . Hence, it follows from ( |4.44| and ( |4.47| that 



(4.47) 



z, 1/z G ri. 



9u 
9n 



(1/r) = -p\^'{1/t)\ = -p 



^'2(t) 



= -P 



W{r)\ 



irr — (r), TGFp. 



dn 



Now, we shall show that (4.481 implies the representation formula 

ds 



u{z) = 



/ log 


Z — T 




1 — ZT 


It, 



\^'{r)\ 



27rp' 



z G flp, 



(4.48) 



(4.49) 



where ds, as before, is the arclength differential. Before we do so, let us gather several 
facts that are easy consequences of the following version of the Green's formula. Let O 
be a bounded open set with smooth boundary dO, and let a and b be two harmonic 
functions in some neighborhood of O. Then 



db , 
a— as 
dO on 



b—ds, 
do on 



(4.50) 
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where d/dn denotes differentiation in the direction of the inner normal of O. In what 
follows, all the partial derivatives are taken with respect to the inner normals of the 
corresponding domains. 

Let L be a smooth Jordan curve and be the bounded compo nent of its comple- 
ment. If a is a harmonic function on some neighborhood of il^ then (4.501 immediately 
implies that 

^"ds = 0. (4.51) 



dn 

Further, denote by the function — log \ z - 



dvz 
dn 



ds 



- Then 

0, z(^nL. 



(4.52) 



Indeed, (4.52) follows from (4.501 applied with a = 1 and b = Vz, by putting O ^ flL in 
the second instance and O ~ CIl \ Dz in the first, where is Dz is any closed disk around 

dvo 



z contained in il^. Thus, (4.501 and (4.52) yield that 



du 

voTrds = 
an, on 



u— — ds = logp 
on, on 



dvp 
dn 



ds 



i/p 



-ds 



'Att log p. 



Moreover, we deduce from (4.48) that 



ri/p 



du 
voTT-ds = 
on 



du 
vo^ds. 
dn 



So, we obtain from (|4.53| and (|4.54| that 

1 



du , 
dn 2 



du 

vo^ds = -27rlogp. 



(4.53) 
(4.54) 

(4.55) 



Finally, let O be fip with the closed disk of radius x around z € flp removed, a = u, 



and b = Vz- Then 

dv 



u— — ds 
do on 



logp 



dvz 



ds — log p 



dn 

27rlogp — 2'Ku(z) 



dvz , 1 

l^ds 

dn X 



uds 



(4.56) 



by (4.52) and the mean- value property of harmonic functions. On the other hand, we get 
du 



that 



VzT:-ds 
do on 



du 
VzT^-ds 
dn 



ri/p 



' dn 
du 

{vz{r)~vz{l/r))—ds 



Vz T^ds — log X 



I du , 

+ / -og^ds. 



log 



1 — ZT 



log 

du 
dn 



du 

I — ZT 



ds 



du 



dn 

ds — 2tt log p, 



ds 



(4.57) 



where we used (4.51), (4.48), and (4.55). Thus, combining (4.56) and (4.57) with (4.50) 
and ( |4.47[ ), we get ( |4.49| . 

Now, we shall construct a triangular scheme such that F is symmetric with respect 
to it. Obviously, this is equivalent to showing that there exists a sequence of sets {-E„}, 
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En = {ej\ri}|"i where D+ is the interior of T, such that |f„| = on F and 

l^nl = 0(1) locally uniformly in Z)+, where 



^«(^) = n f 



Set £1 := const. I^*'!, /p fids = 1. Further, for each n G N, let {F^'"}|"i be a partition of 
Fp into 2n simple arcs that are pairwise disjoint except for possible common endpoints 
and satisfy 



Qds^^, j = l, 

^ o 



Clearly, it holds that 



diam(F^j") < / ds < 



const. 



,2?!. 



,2n, 



(4.58) 



since infr g > 0, where const, is an absolute constant. It also holds that 



k^{T)gds, k^{T) := log 



1 - ZT 



(4.59) 



for any z e F by (4.49). Moreover, for each 2 e F and ti,T2 £ Fp we have that 



|fc.(ri)-fc,(T2)| 



log 



1 



(r2-ri)(l-z2) 



(z - r2)(l - ZTi) 



<0(|r2-Ti|), 



(4.60) 



where the last bound does not depend on z. 

Now, let ej^n be an arbitrary point belonging to F^'", j — 1, . . . , 2n, n E N. Then it 
follows from (4.58 1 and (4.60 1 that for any z G F and t G F^'", j = 1, . . . , 2n, we get 



|fc,(r) - fc,(e,-„)| < 0(diam(F^p'")) < 



const. 



(4.61) 



where const, is, again, a constant independent of z, j, and n. Therefore, we see that the 
functions 

^ 2n 2n „ 



are such that 



< 



(^^(t) - fc2(ej,„)) pds 



E 



const. const. 

gds — 



by (4.59) and (4.61). Finally, observe that |f„| = exp{2nu„} = 0(1) on F. This finishes 
the proof of the theorem since iJ„ C Fp C -D"*" and therefore |f„| = o(l) locally uniformly 
in by the maximum principle and a normal family argument. ■ 
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the following set 



5. Examples 

Denote by i^Q, a e K 

:= 

Clearly, F^ is an analytic arc joining —1 and 1 and therefore is symmetric with respect 
to some triangular scheme by Theorem [T] However, it can be easily computed that 
|r±(e;i)| = 1, t G if e £ F'^ \ {±1}. Hence, if ^? = {E,,}, E„ = {e,,„}2«,, is such 
that Cj^n G \ {^^} fo"^ all n G N and j G {1, . . . , 2n}, then F will be symmetric with 
respect to S'. Further, it can be easily checked that for any e ^ Fa, it holds that 



where 



\r+{e;t)r+{e*;t)\ = \r- {e; t)r- {e* ; t)\ = 1, t e F^, 



2ia 



(1 - q;2) + 2me' 

Thus, if ^ = {^ri} is such that each contains e and e* simultaneously, then F will be 
symmetric with respect to S". Finally, observe also that for a = we get that Fq = [—1,1] 
and e* = e. 

Below, we plot zeros of polynomials orthogonal on F_i/2 with respect to the weights 
Wn{t) = e*t^"(t + 4i/3)^". In other words, these are the poles of the multipoint Fade 
approximants to dfj,{t) = ie*dt/w^{t)\p_^^^, that corresponds to the interpolation 
scheme with half of the points at and another half at — 4i/3. The denominators of 



such approximants are constructed by solving orthogonality relations (2.2 1. Thus, finding 



the denominator of the approximants of degree n amounts to solving system of linear 
equations whose coefficients are obtained from the moments of the measures + 
4i/3)~"(i/Lt(t). The computations were carried with MAPLE 9.5 software using 24 digits 
precision. 




Figure 1. Zeros of qs (disks) and g24 (diamonds). 

In the following example the contour F is generated by ei : 
(87 + 6i)/104, and 63 :~ —i/10, in the sense that 



{i - 3)/4, 62 



{r{ei;t)r{e2;t)r{e3;t))- 



1, 



and is computed numerically. Thus, F is symmetric with respect to any triangular scheme 
such that E^m consists of ei, 62, and 63 appearing m times each, and E^m+i (-£'3^+2) is 
obtained by adding to E^m an (two) arbitrary point (points) from D. Based on the derived 
discretizations of F, the moments of the measures [{t — ei){t — 62) (i — e3)]~^"^h{t)dt, 
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n = 3m, are computed for n = 24 and n = 66, where h{t) = t if lm{t) > and h{t) = i 
otherwise, and the coefficients of the corresponding orthogonal polynomials are found by 
solving the respective linear systems. The computations were carried with MAPLE 9.5 
software using 52 digits precision. 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 



Figure 2. Zeros of 524 (disks) and gee (diamonds). 
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